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Abstract : We present the results of a study performed on the interactions of 10. 6A GeV gold 
nuclei in nuclear emulsions. In a minimum bias sample of 1311 interactions, 5260 helium nuclei 
and 2622 heavy fragments were observed as Au projectile fragments. The experimental data are 
analyzed with particular emphasis of target separation interactions in emulsions and study of crit- 
ical exponents. Multiplicity distributions of the fast-moving projectile fragments are investigated. 
Charged fragment moments, conditional moments as well as two and three - body asymmetries of 
the fast moving projectile particles are determined in terms of the total charge remaining bound in 
the multiply charged projectile fragments. Some differences in the average yields of helium nuclei 
and heavier fragments are observed, which may be attributed to a target effect. However, two and 
three-body asymmetries and conditional moments indicate that the breakup mechanism of the pro- 
jectile seems to be independent of target mass. We looked for evidence of critical point observable 
in finite nuclei by study the resulting charged fragments distributions. We have obtained the values 
for the critical exponents 7, (3 and r and compare our results with those at lower energy experiment 
(1.0A GeV data). The values suggest that a phase transition like behavior, is observed . 



1 Introduction 

The experiment to be described here examines the breakup of relativistic gold nuclei when they 
interact with the target nuclei in nuclear emulsions. The fragments produced are readily identified 
in the emulsions. Specific attention in this paper is directed towards the multiply charged fragments 
that are produced. Some reports on the main characteristics of the interactions have been published 
previously |j|]-|l(|. 

Depending upon the target - projectile combination and the incoming projectile energy, the 
excited piece of nuclear matter decays predominantly by the emission of nucleons, deuterons, tritons, 
helium nuclei and charged particles with 3 < Z < 30 commonly known as intermediate - mass 
fragments (IMF's) and fragments of very heavy charge Z > 31. To understand the dynamics 
involving the formation of helium, IMF's and other multi-fragments in its final state, numerous 
experiments have been performed at low, intermediate and high energies, in both p-nucleus and 
nucleus- nucleus reactions pT|-p^|. 

The experiment to be described here examines the breakup of relativistic gold nuclei when 
they interact with the target nuclei in nuclear emulsions. The relative yields of the different types 
of fragment and the relationships between them are measures of processes that occur during the 
breakup of the excited nuclear remnants. We will study some relationships and compare with those 
observed in others experiments. In this paper we shall present a systematic study on the target and 
projectile fragmentation of the 197 Au - induced emulsion interactions at 10.6^4 GeV from the BNL 
Alternating Gradient Synchrotron (AGS). 

Competing models suggest different decay mechanism and experiments have yet to discriminate 
between several theoretical scenarios which range from the sequential decay of the compound nucleus 
30 , |3T| to statistical nuclear models |32, ^] and percolation models [[34], [35|| . 



It has been stressed out by Campi |36| that the moments of the charge distributions provide a test 
of the suggestion that multifragmentation can be described in terms of percolation theory. If there 
is some critical behavior in the breakup of the nuclei, such as a liquid - gas phase transition, then 
some events should have values of normalized moments much larger than the average. For the low 



energy gold interactions |26| there was a wide range of values and strong and approximately linear 
correlations between the different normalized moments, although it was not possible to conclude 
that then was a phase transition. Also EOS collaboration |27J have reported some of their results 
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from the analysis of 1.0A GeV gold nuclei fragmenting in a carbon target. In their analysis of 
9716 interactions they used the methods developed for determining percolation critical exponents 
to extract the values of specific exponents for nuclear matter from the moments of the fragment 
charge distribution. 

However, one of the main problem encountered in interpreting the results from nuclear emulsion 
experiments is the non-homogeneous composition of the emulsion, which contains both light (H, C, 
N, O) and heavy (Ag, Br) target nuclei. Critical remarks about using minimum biased samples for 
studing critical behavior have been expressed [|16[]. It has been argued that in emulsion experiments 
the mixture between emission sources, both with respect to origin and size, should be a severe 
shortcoming for collisions from few hundreds of MeV/nucleon up to several GeV/nucleon fll6| . 

Therefore one of the main objective of this paper is to present a detailed analysis of specific 
measured quantities for multifragmentation phenomenon for a larger sample of these interactions 
of gold nuclei, with special emphasis on the inclusive interactions with separated light and heavy 
target nuclei. 

Recently EMUOl-Collaboration (JTOj using a statistical analysis based on event by event charge 
distributions, showed that a population of sub-critical, critical and supercritical events was observed 
among peripheral collisions, but the study has been limitted only to critical exponent 7, which in 
some papers p7j is claimed to show little sensitivity to the system under investigation . 

KLMM Collaboration has also looked for evidence of phase changes in the description of 
multifragmentation at 10. 6A GeV, but their results are significantly inconsistent with those reported 
at lower energies, suggesting that percolation theory becomes a less satisfactory representation of 
the breakup for these high energy interactions than it was at lower energies. 

Therefore it will be of great interest to repeat such analysis at the same energy and study the 
values for all critical exponents, using the same methods. 



2 Experimental Details 



Several stacks of electron sensitive, NIKFI BR-2 emulsions, of dimension 10 A 10 X 0.06 cm 3 have 
been exposed to a 10. 6 A GeV 197 Au beam at BNL synchrotron (experiment E863). The sensitivity 
of emulsion was about 30 grains per unit length of 100 fim for singly charged particles with minimal 
ionization. 

Primary interactions were found by along - the track double scanning : fast in the forward 
direction and slow in the backward direction. Fast scanning was made with a velocity excluding 
any discrimination of events in the number of heavily ionizing tracks, slow scanning was carried out 
to find events, if any, with little changed and unbiased projectile nucleus. This analysis immediately 
resulted in a determination of the mean free path (mfp) for interactions. The measured mean free 



path A = 4.99 ± 0.16 cm agreed well with KLMM - collaboration result A = 4.7 ± 0.2 cm p3| and 



an approximation of measured cross-sections on various nuclei and targets A = 4.6 cm . 

In each event the polar angles 9 and azimuthal angles <p were measured. Depending on ionization, 
all tracks emitted from the interaction vertices were classified according to the commonly accepted 
emulsion experiment terminology : 

1) Shower, or s - particles - singly charged particles with a velocity /3 > 0.7 ; they, mainly, 
consist of produced particles and singly charged fragments. 

2) Grey, or g - tracks, whose ionization (the number of grains per unit length ) correspond 
to protons with momentum 0.2 < p < 1 GeV jc ; they consist, mainly of protons knocked 



- out from the target nucleus during the collision with a few percent admixture of it 
mesons with momentum 60 < p < 170 MeV/c. 

3) Black, or b - tracks - mostly, protons with p < 0.2 GeV/c and multiply charged target 
fragments. They have a range of < 3 mm. 

4) Fragments of projectile nucleus - particles with Z > 2. Fragments with Z = 2 are 
identified by visual inspection of tracks, their ionization is constant and equal to g/go ~ 4 
, go - ionization due to relativistic particle (minimal ionization). 

To determine charge of fragments with Z > 3, the density was measured of 5 electrons 
on length no less than 10 mm ; the beam track and track with Z = 2 were chosen as 
tracks for calibrating. The accuracy of charge measurement was for Z < 10, ±1 for 
10 < Z < 28, ±2 for 28 < Z < 40 and ±3 for Z > 40 in units of charge. 

Grey and black tracks amount the group of heavily ionizing tracks Nh = N g + Nj,- 

In each event, the number of produced particles, ir mesons (N n ) have been also deter- 
mined. 

A number of 1311 inelastic interactions were obtained after excluding from the ensemble 
events those of electromagnetic origin and pure elastic scattering. Others details on 
experiment especially on charge fragment measurements and errors, have been recently 
published |IJ. 

3 Projectile breakup 

3.1 Average multiplicity of fast - moving particles. 

The difference between the projectile and the target spect-ator fragments is easy to make. Projectile 
fragments and shower particles are very energetic and they are distributed in a forward narrow 
cone. The angle of this cone is < 6q > — 16.4 mrad and it can be found from the pseudorapidity 
distribution of the shower particles. 

The target fragments are observed as highly ionizing particles, isotropically distributed. They 
can be black particles which are essentially evaporation fragments from the target, with R < 3 mm, 
or gray particles which are knock-out protons or slow mesons with R > 3 mm . 

The breakup of the projectile can be characterized in terms of the numbers, Np - and charges , 
Zp, of the fragments with Z > 3, that are emitted ; the number N a , of alpha particles emitted and 
the numbers N prot , of protons released and available to interact, where from charge balance we can 
write : 

N prot = Z Au - N F * Z F - 2 * N Q (1) 

Z F >3 

We obtain the bound charge by Zb OU nd ,by : 

Zbound = U ( Z ) * Z ( 2 ) 

Z>2 

where n(Z) is the multiplicity of the projectile fragments with Z > 2. 

In Table I , we display results of the present investigation of the average multiplicities < N prot > 
, < N a >, and < Np > of the fragments with Z — 1, Z — 2, and Z > 3, respectively, for the 
197 Au ions. 
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Also a comparison with the results of Cherry et al. [F£J and those of Singh and Jain [[EJ at 
10. 6A GeV are given. A good agreement with both experiments is obtained for < N prot > and 
< N a >, but a slight difference is remarked for mean number of projectile fragments with Z > 3. 

Table 1: The average multiplicity of projectile fragments with Z = 1, Z = 2, and Z > 3 emitted in 
the 197 Au - induced reactions in emulsions at 10. 6A GeV. 



Beam 


Energy 
(AGeV) 


< N prot > 


<N a > 


< N F > 


Ref. 


197 An 


10.6 


28.48 ±0.81 


4.63 ±0.13 


2.01 ±0.06 


Hi 


197 Au 


10.6 




4.53 ±0.13 


1.91 ±0.06 


Hi 


197 Au 


10.6 


28.44 ±0.64 


4.51 ±0.08 


2.37 ±0.03 


This work 



The numbers of fragments emitted from individual interactions can be compared by looking at 
the distributions of the numbers of events with a given number of fragments divided by the total 
number of events, the fractional yield. These fractional yields for fragments, those with Z > 3 (Fig. 
la), for alpha particles (Fig. lb), for heaviest fragments (Fig. lc) and for numbers of released 
protons (Fig. Id) are compared in Figure 1 with the published results of KLMM collaboration 



23] , p4| . It appears that these yields are nearly identical, although some differences appear especially 
for heaviest fragment and number of proton released. 

The experimental spectra for number of released protons (Fig. 2a) and number of produced 
pions (Fig. 2b) as well as their correlation with number of mesons (Fig. 2c) and the number of 
heavily ionizing tracks (Fig. 2d) are given in Figure 2. As would be expected, there is a strong 
correlation between the number of proton released and the number of pions produced - there is a 
rather well defined upper bound to the number of pions N n for a given number of proton released 
Nprot. Fig. 2d shows that interactions with light target nuclei, those with N h < 7, never results 
in a copious pion production. For Nh > 8 , where the target must have been a AgBr nucleus, N n 
appears to be almost independent of Nh- 

It is interesting to note that our data set for 197 Au do not give any evidence of the occurrence 
of binary fission in the charge range of 35 < Z < 45 |§ , although a significant enhancement of the 
fission events was observed at lower energies with the 197 Au projectile, as reported in Ref. |2"6f. 



3.2 Correlation between N a and Np 

Broad characteristics of the projectile fragmentation can be explored through the correlation be- 
tween the multiplicity of helium particles Ac, and heavier fragments Np . Such a correlation is 
given in Table 2 for the 197 Au data . From Table 2 one can see that the number of pure central 
interactions, in which the projectile has disintegrated completely into singly charged fragments 
with no emission of heavier fragments, is very small (~ 1 %). More than 20 % of the events have 
N a = and Np = 1. Only 2.5 % of the total number of interactions is found with Np > 5 and 
< N a < 10. On the other hand, only 2.1 % of the total umber of collisions is observed with 
Np < 5 and N a > 10. No event is observed with Np > 8 and N a > 10. The maximum number of 
fragments with Z > 3 in our sample is Np = 8 and that of fragments with Z > 2 is N a ± Np = 21. 



Table 2: Characteristics of projectile fragmentation of the 197 Au beam through a correlation between 
the number of helium particles, N a and the number of heavy fragments, Np with Z > 3. 
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Specific correlations for mean number < Np + N a > with number of singly charged fragments 
(Np) are shown in Fig. 3a. We can remark an expected linear correlation for Np < 10 and a plateau 
for Np > 10. Correlations with number of shower particles N$ are given in Fig. 3b, ( we note that 
N s does not include singly charged fragments Np) indicate a clear change in slope for N s > 70, 
which is due to separation of central and peripheral collisions. A rise and fall for correlations of mean 
number < Np + N a > with number of released proton, which are complementary to correlations 
with Zfound (see also Fig. 5a,d) are depicted in Fig. 3c. A more complex correlations are obtained 
by plotting < Np + N a > as function of number of black and grey tracks (Nh) in Fig. 3d, in which 
we can see a different slopes for interactions with light (N h < 7) and heavy nuclei (Nh > 8. 

3.3 Target separation 

Nuclear emulsion are a composite medium composed of AgBr, CNO and H. Certainly, they are also 
other nuclei in emulsions, but their concentration are too small to be taken into account. It was 
a difficult task to separate interactions on different classes of targets and it is impossible to find 
certain separation criteria that give no admixture between those classes although there are many 
correlations between the measurable parameters that give informations regarding the target nuclei. 

The separation technique we used was based on the analysis of specific correlations between 
target break-up and particle production. Depending upon the target break-up, we separate the 
sample of 1311 interactions into three classes : 

a) N h — 0, 1 ; this class includes all Au + H interactions but also interactions with other targets. 



b) 2 < Nh < 7, containing the rest of Au + CNO interactions not included at a), but also some 
Au + AgBr. 

c) Nh > 7, including only interactions with AgBr. 

An additional relation that give some target separation for a) events is based on the distribution 
of the number of shower particles with 9 > 9 . Events with no black or gray tracks include most of 
the H interactions, but also the most central events with CNO and AgBr. The distribution of the 
shower particles for these events indicates the limit of no more than about 40 relativistic particles 
with 9 > 9q emitted from Au + CNO interactions. The separation between H and CNO peaks give 
an admixture of less than 10 % for this class of events. 

For Nh = 1 we analyzed separately elastic interactions from a kinematic approach. Knowing 
the emission angles for the black track and the recoiled gold nucleus we could compute the range- 
momentum ratio for most of the elastic events. This ratio is highly dependent on the target size 
and it's very useful in order to separate elastic interactions with CNO from those with AgBr. For 
elastic interactions with hydrogen nuclei, the recoil proton appear most often as a gray track. The 
number of shower particles for inelastic interactions with Nh = 1 indicates that AgBr events are for 
N 8 , > 40. 

Class b) events were separated by plotting the number of shower particles against N^. The 
separation between CNO from AgBr populations is given by : 

N a > < 175 - 14.5 • N h (3) 

where, N s > = N s + N P . 

The reason we took the limit Nh > 7 for class c) events was that 8 target tracks for an oxygen 
nucleus correspond to a total charge break-up, reasonable to assume that at least one of the released 
target protons become a relativistic particle. 

We plotted the distributions of the number of a particles and PF's for the interactions of gold 
nuclei with HCNO and AgBr targets in Figure 4. We can see in Fig. 4a that the emission of 
a particles is enhanced for interactions with heavy targets, the distribution being shifted to the 
right. However the distributions of heavy fragments in Fig. 4b show the same behavior for the two 
samples. 

The correlations between the mean number of emitted alpha's and fragments and the total bound 
charge are depicted in Fig. 4c and Fig. 4d. We can see that for peripheral events corresponding to 
Z bound > 50 the emission of alpha's and fragments have the same behavior both for light and heavy 
emulsion targets. However a particle emission is suppressed almost down to zero when centrality 
is increasing, especially for HCNO interactions. The difference between the two samples for PF's 
is not so obvious. Nevertheless we can notice some decrease of the mean number of fragments for 
HCNO compared to AgBr in semi-central events. 



4 Specific Correlations for Multifragmentation 

4.1 Correlation between some observable and Zh OU nd- 



The parameter Zb oun d, which was defined by Hubele et al., |19| is related to the size of the projectile 
spectator and the energy deposited in a given collision can also be explored through Z bound . Z bound 
is always larger than or equal to Z max . 

The heaviest fragment charge in each collision, Z max , provides useful information on the exit 
channel of that collision. Fig. 5a shows a scatter plot of the correlation between Z max and Z boun d 
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for individual events of the 197 Au data. We remark that most of the data points are situated below 
the diagonal, that means that Z bound is always larger than or equal to Z max . Also in the peripheral 
collisions, the largest fragment contains most of the total bound charge and in such reactions Z max 
may be identified as the heavy residue of the beam nucleus after evaporation. For the central events, 
we observe that Z max becomes a smaller fraction of Z^^ . This plot also confirms that symmetric 
fission is a very rare kind of events. 

In Fig. 5b we give the variation of mean value of < Z max > /Zb ea m as a function of Zbound for 
the 197 Au data. A sharp rise of this ratio can be seen from Fig. 5b for Zbound ~ 40 — 79. For 
Zbound ~ 2 — 40 the rise is rather slow. This plot give some information into the degree of breakup 
of the projectile nucleus. 

In Fig. 5c we plot the average multiplicity distributions of alpha particles < N a > and in Fig. 
5d the mean number of intermediate - mass fragments (IMF's) with 3 < Z < 30 as a function 
of Zbound for 197 Au projectiles. Our experimental data confirm, the so called "rise and fall" of the 
multi-fragment emission . 



4.2 Charged fragment asymmetries 

In this chapter, we investigate the charge asymmetry between the largest fragments (Z max ) and the 
second largest Z F M2 fragment in an event, using the two body relative asymmetry Ras, defined as 



19 



„ ^max ~ ^FM2 , . n. 

Ras = ~ —7} (4) 

£>max T £>FM2 

Only fragments with Z > 2 are included. 

Also Rasi, the asymmetry between the second largest Zfm2 and third largest Z F M3 in an event 
is obtained by : 

Z FM2 — Z FM3 

Rasi = -= — ^ (5) 

^FM2 + ^FM3 

We can also investigate the breakup process in a more qualitative manner through another 
parameter known as three body asymmetry [19[] as : 



'{DFAf + {DFBf + {DFCf 

Ras2 = wtj^ (6) 



where : 



DFA = 


Zmax 


<z> 


DFB = 


ZfM2~ 


< z > 


DFC = 


ZfM3~ 


< z > 



with : 

< Z >= -{Z max + Z F M2 + Z FM3 ) (7) 

The parameter Ras2 has a maximum near unity when there is heavy residue of the projectile 
spectator and has a zero value when the projectile fragments of equal size are emitted in the 
collisions. 
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For the analysis of two body asymmetries given by Eq. (|]) and Eq. (^) , we have selected the 
events with at least three fragments of Z > 2. 

In Fig. 6a and Fig. 6b, we plot a correlation between mean value of ratio < Ras > an d 

< Rasi > respectively with Z bound for 197 Au data. The ratio < Ras > decrease monotonically from 
its maximum value ~ 0.9 to almost zero as one approaches from extremely peripheral toward more 
violent collisions. The ratio < Rasi > ^ises linearly with Z bound for Z bound < 40 and decreases for 
40 < Zbound < 79 for these class of events selected (Np + N a > 3) . In Fig. 6c we represent ratio 

< Ras > for IMF's and we see that the ratio remain almost constant for 20 < Zbound < 60. For 
Zbound > 60 a clear leading fragments effects appear. In Fig. 6d we plot a graph of a correlation 
between mean value < Ras2 > and Zb oun d- The parameter < Ras2 > rises almost linearly with 
enhancing value of Zb oun d- 

In order to observe if there is any difference in the behavior of the projectile break-up mechanisms 
for interactions with different targets in emulsion, we analyzed two and three-body asymmetries 
between the heaviest emitted fragments. 

In this analysis we used together the separation criteria for H and CNO interactions in order 
to eliminate any possible admixture between these groups. Thus we built two samples, the first 
containing interactions with light (HCNO) nuclei and the second AgBr events. 

In Figure 7 we plot the same parameters for interactions with light target (HCNO) (open circles) 
and heavy targets (Ag,Br) (full circles) as described in section 3.3. 

We plotted the correlations between the mean two and three-body asymmetry ratios for inter- 
actions with light and heavy nuclei and the total bound charge. We can see in Fig. 7a that the 
asymmetry ratio between the first and second heavy fragments increases monotonically both for 
light and heavy samples from most central interactions to peripheral one [Zbound — 79). The 
second and third fragment asymmetry ratio (Fig . 7b) seems to be constant for semi-central and 
peripheral interactions, decreasing down to for Zbound < 40 when increasing the centrality. 

The asymmetry ratio for IMF's (see Fig. 7c) is constant in a wide region of Zbound, peaking near 
the value of 70 which indicate a higher asymmetry between the IMF's in peripheral interactions. 
The three-body asymmetry ratio plotted in Fig. 7d show an approximately linear dependence with 
Zbound but also there is no significant difference between the two samples. 

These observations point out that multifragmentation of gold nuclei at this energy does not 
depend on the target nucleus. Thus it doesn't matter that nuclear emulsions contain an admixture 
of targets from the multifragmentation mechanisms point of view. 

4.3 Moments of the Charge Distribution 

Multifragmentation has been considered to be one of the most important aspects of heavy - ion 
collisions since it has been speculated that the decay of a highly excited nuclear system might carry 
information about the equation of state and the liquid-gas phase transition of low density nuclear 
matter. The similarity between statistical multifragmentation models and percolation theory has 
been stressed in ref. P4"| , |33| The relevance of percolation ideas in nuclear fragmentation can be 
investigated better by examining cross relations between various moments of the fragment size 
distribution. We will show in this paper that experimental data have strong similarities with the 
predictions of percolation models. 

Following Campi suggestion we investigate the moments of the charge distribution of the nz 
projectile fragments (PF's) using an event-by-event based analysis. For each event, we determine 
the multiplicity of charged fragments, nipp, and the number of charged fragments, nz, of nuclear 
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charge Z. We then consider the i'th moments of this distribution 



M i = Y t Z i *nz (8) 
z 

where the sum is extended over all the fragments except the biggest cluster (fragment) which is 
being considered as the percolating cluster f39f . Physically Z max corresponds to the bulk liquid in 



an infinite system. 

The zero order moment is obtaining by taking % = in Eq. (P) : 



M = J2n z (9) 



It has also been suggested [36| that the conditional moment, 72, which is combined from the 
moments M 2 , M and Mi as : 

M 2 M 

72 = -mF (10) 

give more selective information. 

In Figure 8 the mean values of M (Fig. 8a), M\ (Fig. 8b), and M 2 (Fig. 8d) , averaged 
over small range of Z^und for the nonzero values of the moments, are shown as function of Zbound 
. In Fig. 8c we represent a variation of mean values of < 72 > for the events with at least three 
charged projectile fragments with Z > 2. The value of < 72 > increases rather slowly in the range 
of 2 < Z bound < 50. The maximum value of < 72 > ~ 1.35 at Z bound ps 50 and then decreases 
for Zbound > 50. Some fluctuations in < 72 > can be seen in this region for the 197 Au data. For 
the values of Z bound < 20, < 72 > = 1 (events in which the projectile fragments of the same 
size are emitted). For the infinite nuclear system, the scaling theory of critical phenomena predicts 
that at the critical point 72 diverges at a rate that depends on the critical indices of the phase 
transition [ 36| . For a finite nuclear system, 72 is predicted to have a smooth behavior ||36|| . From 



Fig. 9c we remark a smooth increase of < 72 > up to Z^und ~ 60 for the 197 Au data. However 
some fluctuations appear in our 197 Au sample . Similar results were analyzed for 197 Au data at 
10.6A GeV and 208 Pb data at 160A GeV |Z§. The peak observed is due to finite size of the nuclear 
systems under investigation. 

In Figure 9 we study the dependence of the same quantities M (Fig. 9a), Mi (Fig. 9b), M 2 (Fig. 
9d) and 72 (Fig. 9c) function of Z^und for different interactions on light targets (open circles) and 
heavy targets (full circles). No significant differences between the two target groups are observed. 



5 Critical Exponents in Multifragmentation 

5.1 Leading fragments 

Fig. 10a, exhibits an area plot of the correlation between the charge of the heaviest fragment Z max 
and the first order moment M\ ( or remaining bound charge Z r ) . This shows that there are many 
interactions where Mi > Z max and that, even when there is a well defined leading fragment (a 
fragment that carries more than half of the bound charge Zb OU nd) - it is usually accompanied by 
appreciable other bound charge. 

Zbound Z max -\- M^i Z max -\- Z r (H) 

1 n 



The plot of Z max versus M±, already presented in Fig. 10a shows a clear, but complex correlation, 
which suggest two distinct populations. This suggestion is enhanced by separating the interactions 
into those with ( M x < Z max ) and without (Mi > Z max ) leading fragment. 

If the mean values of < Z max > are determined over small intervals of Mi, for these two classes 
they will appear well separated (Fig. 10c). 

A similar comparison can be drawn from the correlation between the next heaviest fragment 
Zfm2 and Z max in Fig. 10b. The mathematical restrictions are indicated by the dasheded lines on 
this plot, but it is clear that the allowed space is not uniformly populated. In particular the fission 
region near the apex is almost empty (see also reference f9|), in contrast to the situation seen at 
lower energies ( [f26fl ). 

Considering the ratio R\ for all fragments, where : 

D Z max 

1 = ^T (12) 

Then there is a leading fragment, as defined above , when R\ > 1 . We see from Fig. lOd that 
ratio Ri is a strong function of Zb OU nd- Thus, while there is a well defined leading fragment when 
Z bound is greater than 60, it is much less well defined when Z bound falls to 40. 

5.2 Critical Exponents 

The observation of a power law behavior for the size distribution of the fragments has triggered 
a number of studies that have looked for evidence of critical behavior pTfl , pjfl , p2| , flQfl . These 



analysis consider nuclear multifragmentation as one example of a critical phenomenon and attempts 
are made to extract from the data the related critical exponents. 



The EOS Collaboration |27] have reported some of their results from the analysis of 1.0A GeV 
gold nuclei fragmenting in a carbon target using the methods developed for determining percolation 
critical exponents to extract critical exponents for nuclear matter from the moments of the fragment 
charge distribution [^] : 

We assume that the multiplicity of fragments, m = Np + iV^ + N prot is a linear measure of the 



distance from the critical point as suggested by Campi pM. Here Np - is multiplicity of fragments 
[Z > 3, Na is multiplicity of alpha particles [Z = 2 and N prot is the number of released protons 

{Nprot — Zbeam ~ Zbound with Zbeam = 79). 

The region in m below the assumed critical multiplicity m c is designated as the liquid phase and 
that above m c as the gas phase. It is assumed that in the liquid phase the heaviest fragment Z max 



is omitted in forming the moments, but is not omitted when in the gas phase |39[ . 

Also our analysis tacitly assumes that all of the projectile-related charges are associated with 
multifragmentation and include in the analysis of the moments number of released protons, N prot 
as fragments, modifying Eq. (§) and replacing Mj by M*. 

The critical exponents 7, /3, and r for large systems are given by the following equations in 
terms of the multiplicity difference, ( = m — m c by : 

Ml ~ |C 7 | (13) 
Z max ~ \C\ P (14) 



n z ~ Z- T (15) 
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These exponents 7, (3 and r are related by the equation 



- 2 + jh < 16) 

We start our analysis by examining the variation of the mean values for second moments < > 
as a function of multiplicity m which is depicted in Fig. 11a . A relatively abrupt change in this 
distribution is apparent for m c around 26, suggesting that there could be a phase change at this 



critical value of m c . This values is similar to those reported by others experiments |27j , f23[ . 

Fig. lib shows a log-log plot of mean values of second moments < M| > with ( for the assumed 
liquid and gas phases, setting m c = 26 ( where log stand for natural logaritms). The clear separation 
between the two phases arises from the inclusion of Z max in the determination of < M| > in the 
gas phase. In Fig. 11c and Fig. lid the values are represented in a scatter plot. 

If we examine over the entire available range of ( , neither phase shows the power law behavior 
predicted by Eq. [13] . However, if rather narrow regions of \(\ are selected, 5 < \(\ < 20, then 
a good fit to such a power low can be obtained for the gas phase , with 7 flas = 0.86 ± 0.05 and a 
reduced x 2 °f 2.94. A fair fit can also be obtained for the liquid phase, with "/liquid — 0.83 ± 0.14, 
reduced x 2 °f 1-44. The results are represented in Fig. 12a. 

We note that the values are relatively sensitive to the range of |£| used, as we expected, since 
finite size distorsions dominate as \(\ — > 0, and signature of critical behavior vanish for large 
corresponding to mean field regime. Adding two more values to those shown in Fig. 12a changes 7 
to = 0.69 ± 0.11 with x 2 of 1.56 and ^ gas = 0.73 ± 0.04 with X 2 of 4.3. 

The most important results is that the values for juq U id are close to the values of 7 9as , which 
implies that the conditions for a phase change have been satisfied. No better match can be found 
for m c = 30. In practice to estimate the uncertainties, we varied the fitting region by changing the 
upper and lower multiplicity limits. The overall estimated uncertainties are 14.5 %. 

We continue the analysis determining the exponent (3 from Eq. [14| considering the liquid phase, 
where Z max is well defined. Fig. 12b shows log(< Z max >) as a function of |£| . The value obtained 
for (3 = 0.25±0.02 with a reduced x 2 = 1-56 is in a good agreement with that of 0.29 ±0.02 reported 
at 1.0 A GeV |27j (see also table The value obtained for an exponent is sensitive to the range of 
values chosen for \(\ . 



The critical exponent r in Eq. (15) can be determined from the slope of log(< M| >) versus 
log(< M% >) (see Fig. 12c) using only the gas phase Fig. 12d shows a power law fit with 
r = 2.23 ± 0.05 and a reduced x 2 of 1.51, for 3 < |C| < 36. 

1 1 or 



Our value is practically the same with the value predicted in infinite percolation models [35 



the values reported for the 1.0A GeV |27| 2.26 or the value of 2.23 calculated from our measured 
values of (3 and 7 gas using Eq. ([16]). 

Reducing the range of \C,\ values used for this fit does not make any significant change in the 
deduced values of r. We note that the critical exponent r is close to 2.2 for many three-dimensional 
systems and thus does not permit a determination of the universality class of phase transition. 

The exponent values are summarized in table |3|, which also list the values from lower energies 
experiment [p7H , percolation [j39fl , liquid-gas values and the mean field limit of the liquid-gas 
system [f|I . 
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Table 3: Critical multiplicity and exponents for Au projectile fragmentation and other 
three-dimensional systems 



Quantity 


m c 


7 


P 


T 


Our Exper. 


26 


0.86 ±0.05 


0.25 ±0.02 


2.23 ±0.05 


EOS Exper. j27|| 


26 ±1 


1.40 ±0.1 


0.29 ±0.02 


2.14 ±0.06 


Percolation |39| 




1.8 


0.41 


2.18 


Liquid - gas 
mean field 41] 




1.0 


0.50 


2.33 


Liquid-gas |40|| 




1.23 


0.33 


2.21 



We note that the values of 7, j3 and r obtaining using this method obey the scaling relation, 
Eq. flip]). By varying the fitted region, we have obtained exponents which differ by as much as 
15 %. We can take this to be a measure of the uncertainty in their values. It is significant that the 
values for (3 and 7 are different from either the percolation or the mean field values. Also we can 
remark that 7 exponent is different from those of nominal fluids and depends on energy. 

6 Conclusions 

In the present work we have studied, the properties of the projectile associated particles emitted in 
interactions of the 197 Au ions accelerated at an energy of 10. 6A GeV obtained from the BNL AGS. 

The average multiplicities of the fast - moving projectile particles such as < N a >, < Np >, < 
Nprot > seem to depend upon the mass of the target. The majority of the multiply charged fragments 
are helium nuclei, while the majority of those fragments with Z > 3 are light. The multifragment 
emission is a dominant reaction channel as observed when the distributions of < N a >, < N IMF > 
are represented as function of Z bound . These distributions are peaked at Z bound « 35 — 40 and 
shows slight dependence on the target mass. 

Nuclear emulsion detectors provide an excellent tool to study the global characteristics of 
nucleus-nucleus interactions since they allow a simultaneous investigation of the processes of nuclear 
fragmentation and multiple particle production and allow a study of the correlations between these 
processes. 

Even if the emulsion detector contains different targets, multifragmentation, when expressed in 
terms of Z bound , appears to be relatively insensitive to the nature of the target and the results can 
be compared with those from studies using pure targets. 

An analysis of the moments M Q , Mi and M 2 as well as conditional moments such as 72 also 
proves that the breakup mechanisms has no dependence on the target size and a broad peak in the 
72 - Z bound relation shows that the nuclear systems employed in the present investigation induced 
a finite size effect. 

Two and three body asymmetries are explored through the distributions of < Ras >,< Rasi > 
and < Ras2 > as a function of Z bound . Within statistical errors, the distributions shows almost 
similar behavior on different target which indicates that the breakup mechanism has no dependence 
on target mass. 

A study of multiplicities suggest that there could be a phase change in the residual nucleus 
that depends on the multiplicity of the charge fragments, in a manner similar to that predicted by 
theories such as percolation that describe the process of multifragmentation. 

1 3 



Our analysis for a critical point and a phase change based on our high energy data give results 
which are consistent with the analysis reported for the low energy results. The presence of a critical 
point is well established from our data. Comparison with percolative and liquid-gas systems show 
remarkable similarities. However, some essential differences on values of critical exponents are 
pointed out. 

To further characterize this phenomenon, we must determine whether the fragmenting system 
is thermalized and if so its temperature and density [Q, and also whether the multiplicity is 



proportional to temperature. Some results have been published recently J43], [44 
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Figure 1: Part (a) - Fractional yields of fragments with Z > 3 observed in the two experiments ; 
EMU-01 data (full circles) and KLMM data (open squares). Part (b) - Fractional yields of alpha 
particles in the two experiments as a function of alpha particles emitted from an interaction. Part 
(c) - Fractional yields of the heaviest fragments in the same experiments as a function of the charge 
of the fragment. Part (d) - Fractional yields of released protons as a function of the numbers of 
protons released from an interaction. 
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Figure 2: Part (a) - Distribution of number of proton released - full circles EMU-01 data. Part (b) - 
a comparison of the number of pion produced in the two experiments. EMU -01 data (full circles) ; 
KLMM data (open squares). Part (c) - Correlation between the number of protons released and the 
number of pion produced . Data from EMU-01 experiment. Part (d) - Number of pions produced 
as a function of the number of black and grey tracks (N h ) emitted from the target. 
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Figure 3: Correlation of mean number of alpha particles and fragments with Z > 3 ( < N a + Np >) 
with number of singly charged fragments (iVp) (Fig. 3a), with the number of shower particles (Ns) 
(Fig. 3b), with the number of released protons (N prot )(Fig. 3c) and with the number of black and 
grey tracks (N h )(Fig. 3d). 
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Figure 4: Part (a) - The distributions of the numbers of alpha particles in interactions with a 
emulsion light (HCNO) ( open circles) and heavy (Ag,Br) (full circles)- EMU-01 data. Part (b) - 
The distributions of the numbers of fragments with Z > 3. Part (c) - The average number of alpha 
particles < N a > as a function of Z bound . Part (d) - The average numbers of fragments ( Z > 3 ) 
as a function of Zb oun d- The experimental data have the same meanings as in Part (a). 
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Figure 5: Part (a) - A scatter plot between the correlation of largest charge Z max and Z boun ^ 
in individual events . The diagonal is shown by a dashed line. Part (b) - Mean values ration 

< Z max / Z beam > as a function of < Z bound >. Part (c) - Mean numbers < N a > as a function of 

< Z bound > for full sample. Part (d) - Mean numbers of intermediate fragments (IMF's) < N IMF > 
(with 3 < Z < 30) as a function of Z boun d- 
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Figure 6: Part (a) - Two body relative asymmetry < R^s > versus Abound- Part (b)- Two body 
relative asymmetry < Rasi > versus Zb OU nd- Part (c) - Two body relative asymmetry for IMF's 
versus Z bound . Part (d) - Three body asymmetry < Ras2 > versus Z bound . For the definitions of 
charged particle asymmetries < Ras >, < Rasi >, < Ras2 > see section 4.2 
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Figure 7: Part (a) - Two body relative asymmetry < Ras > versus Zf, ound . Part (b)- Two body 
relative asymmetry < Rasi > versus Zbound- Part (c) - Two body relative asymmetry for IMF's 
versus Z bound . Part (d) - Three body asymmetry < Ras2 > versus Z bound . Experimental points have 
the same meanings as in Figure 4. For the definitions of charged particle asymmetries < Ras > , 
< Rasi > , < Ras2 > see section 4.2. 
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Figure 8: Mean values of zeroth moments, < M > (Fig. 8a), first moments, < Mi > (Fig. 8b) 
and second moments < M2 > (Fig. 8d) as a function of Zb OU nd- Part (c) - Variation of the mean 
values of the conditional moment < 72 > with Z bound . For the definitions see section 4.3 
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Figure 9: Mean values of zeroth moments, < M > (Fig. 9a), first moments, < Mi > (Fig. 9b) 
and second moments < M2 > (Fig. 9d) as a function of Zb oun d. Part (c)- Variation of the mean 
values of the conditional moment < 72 > with Z bound . The experimental data points have the same 
meanings as in Figure 4. For the definitions of the symbols see section 4.3 
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Figure 10: Part (a) - Correlation between charge on heaviest fragment Z max and remaining bound 
charge Z r = Mi. Right diagonal dashed line show charge limit, left full line shows where Mi = Z max . 
Part (b) - Correlation between heaviest, Z max , and second heaviest fragments. Diagonal dashed 
lines show charge limits, left full line shows where Z FM2 = Z max . Part (c) - Mean of the charge 
on the heaviest fragment < Z max >, as a function of the first moments, Mi separated into those 
events where Z max is greater or less than Mi. Part (d) - Mean values of the leading fragment ratio 
< Ri > = Z max /Mi as a function of bound charge 
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Figure 11: Part (a) - Mean second moments < Mjf > as a function of multiplicity m. Part (b) - 
Mean second moments < M| > as a function of the multiplicity difference ( assuming m c = 26. 
Part (c) - A scatter plot of second moments M| as a function of of the multiplicity difference ( for 
liquid phase. Part (d) - A scatter plot of second moments as a function of of the multiplicity 
difference |C| for gas phase. 
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Figure 12: Part (a) - A log-log plot of mean second moments < M| > as a function of multiplicity 
difference ( assuming m c = 26, over a limited range of (. Linear fits are shown for both the gas and 
liquid phases to define the exponent 7 (see Eq. |13|). Part (b) - A log- log plot of the mean values of 
heaviest charge, < Z max > , as a function of multiplicity difference (, in the liquid phase, assuming 
m c = 26. A linear fit is shown to define the exponent (3, Eq. [14|. Part (c) - Correlation between 
mean second and third moments < M| > and < M| >, for the gas phase, assuming m c = 26. Part 
(d) - A log - log plot for correlation between mean second and third moments < M| > and < M3 >, 
for the gas phase, assuming m c = 26 over a limited range of (. A linear fit is shown to define the 
exponent r, Eq. [L5|. 



